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Abstract 

Numerous models have been proposed to describe the Bose-Einstein 
correlations in multiple particle production process. In the present pa- 
per we describe a generalization, which includes many previous models 
as special cases and, therefore, can be useful for work of comparison. 
We apply the powerful methods of eigenfunction expansions and gen- 
erating functionals, which often make the calculations much shorter 
than in the original papers. 



1 Introduction 

Bose-Einstein correlations in multiple particle production processes at high 
energy are now much discussed. The reviews JI], ||, || contain hundreds of 
references and many more can be found in the papers quoted there. There 
is a variety of models often based on very different physical pictures of the 
production process. In the present paper we stress that many of these models 
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differ only in their choice of a single particle density matrix, further called 
input single particle density matrix. Once this choice has been made, the 
calculations become model independent. Following the approach from |^, |^ 
we describe a general model with a free function, which contains all these 
models as special cases. Then we show that the further calculations can be 
greatly simplified, if one uses eigenfunction expansions for the input single 
particle density matrices and suitable generating functionals to summarize 
the information about the distributions, which can be compared with experi- 
ment. This paper is essentially a description from the point of view presented 
in 1^, Q of a class of models including all the multiparticle symmetrization 
effects. 

2 Models with factorization 

Many models of Bose-Einstein correlations in multiple particle production 
processes can be reduced to the following four steps, which we list here post- 
poning their discussion to the following sections. 

1. Guess an input single particle density matrix Pi'\p,p'). 

2. Construct the n-particle density matrix for distinguishable particles 

n 

Pn{PU--^Pn;P'l, - ■ ■ ,P'n) = H Pl\Pk, P' k) ■ (l) 

k=l 

3. Symmetrize this density matrix in order to obtain a density matrix for 
undistinguishable particles. The diagonal elements, which are enough 
to calculate all the momentum distributions, are 

PniPl, ■ ■ ■ ,Pn) = ^Y.Pn\P<rl, ■ ■ ■ ,Pan]PTl, ■ ■ ■ , Prn) , (2) 

where each of the summations over a and r extends over all the n\ 
permutations of the set of indices 1, . . . , ra. The normalization is such 
that calculating the trace one integrates over the momentum space 
without introducing the factor —.. 

4. Build the diagonal elements of the overall density matrix according to 
the formula 
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oo 



p = ^pW(n)p„, (3) 



n=l 



The coefficients are usually chosen Poissonian with some average v. 

n 

p(o)(n) = ^e-^ (4) 



Further, models built according to this recipe will be called models with 
factorization. Let us now discuss the four steps. 

3 Step 1 — choice of the input single particle 
density matrix 

The superscript of the function pf^ indicates that this is not the observed 
one-particle density matrix, but an input function necessary for the construc- 
tion of the true density matrix. It may be interpreted as the single particle 
density matrix for the unphysical case, when the particles are produced in- 
dependently and there is no Bose-Einstein symmetrization. The different 
models differ in the inspirations used to guess this function. The shortest 
way is, of course, to guess directly the input density matrix pi \ however, a 
longer way making use of a model may be easier. Let us quote some examples. 

One can guess a source function S{X, K) and calculate the density matrix 
pf^ from the formula 

pf\pj) = I d^Xe-^^^'SiX, K), (5) 

where 



K = ^{p + p') (6) 
q = p-p' (7) 

are four- vectors with the four- vectors p, p' being on shell particle momenta. 
The advantage of this approach is that the source function reflects the space- 
time and momentum distribution of the sources of particles. Thus, there are 
intuitions what it should look like 0]. 
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Another approach is to assume that the particles originate from a large 
number N of independent, incoherent sources. In order to produce a reason- 
able number of particles, each of the sources must be weak. A formula for 
the "total source" p, Q with a good high limit is 

1 ^ 

J(p, AT, C, 0) = ^ E e*^'=e^^^'=jo(A.p), (8) 
vA/ 

where c/) = {</)i, . . . , ip^} is a set of random phases, ( = {xi, vi, . . . , xjy, v^} is 
a set of parameters characterizing the A^ sources and A^ are Lorentz transfor- 
mations with velocities Vk, acting on the momentum fourvector p. The "total 
source" is not a source in the sense of quantum field theory. It is rather a 
kind of single-particle wave function with the condition that the components 
labelled by different indices k are not allowed to interfere. This is imple- 
mented by the random phases (pk- An object of this kind can be replaced by 
a single particle density matrix 

1 ^ 

p'Mp', C) = T7 E e''-^^-^yo{AkP)fo{AkP')- (9) 
fc=i 

There are two ways of proceeding further. One can assume that the density 
matrix is an average of this density matrix over the sets (k = {xk, ^4} and 
obtain 

pf\p,p') = J dCkp{Ck)e''''^''~''^Jo{AkPmAkP'). (10) 

In the approach of ref. ||^ there are additional factors dependent on the 
number of sources A^, which for A^ large drop out from the final result. We 
absorb the A^-dependent factor into the density p(Cfc) so that it does not 
appear explicitly in the formulae. In this approach one has to guess the 
spectrum of a single source at rest jo{p) and the distribution of sources p(Cfc)- 
The resulting single particle input density matrix does not depend on the 
unmeasurable parameters any more. 

Alternatively one can postpone the averaging over the unmeasurable pa- 
rameters as done in ref. (further quoted GKW) and choose the pure state 
input single particle density matrix 

p^'^iP'P') = .^ Jip,N,cA)riP',N,CA), (11) 

n A, C,<P 
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where 



n{N,CA)= d'k\J{k;N,CA)\' 



(12) 

is a normalizing constant. It is natural to assume in such models that also 
the input multiplicity distribution (n) depends on the unmeasurable pa- 
rameters. Models of this kind do not belong to the class of models discussed 
in the present paper, unless simplifying assumptions are made. We will dis- 
cuss only the simplest case, when the unmeasurable parameters are fixed or 
absent, and consequently no averaging is necessary. This model (cf. [Q) will 
be called pure state model, because for it the input single particle density 
matrix, which can be written in the form 



(13) 



corresponds to a pure state. There are two reasons to consider this grossly 
oversimplified model: it is one of the very few models, where multiparticle 
effects can be included analytically and it is a good starting point for the 
discussion of the much more important GKW model. 

Still another strategy is to guess a set of single particle wave packets \a) 
and the distribution of such packets p{a) |T^, |TT|. Then 



Pi\p,P') 



da{'p\a) p{a) {a\p') . 



(14) 



Let us note for further reference that if the states \a) form an orthonormal set, 
they are the eigenvectors corresponding to the eigenvalues p{a) of the matrix 
pT'iPiP')- Then, usually, the integration gets replaced by a summation. 

Whatever the starting point, very often one finally obtains a Gaussian 
(cf. e.g. i i 0) 



Pi\p,P') 



n 

i=x,y,z 



27rA2 



: exp 



2A 



(15) 



In order not to contradict the Heisenberg uncertainty principle, one must 
have 



(16) 
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In the wave packet approach the parameters i?j and Aj are expressed in 
terms of other parameters in such a way that this condition is automatically 
fulfilled. In other approaches one must impose it as a constraint. 

One could include in the exponent a term isq ■ K [jl2| , where s is a real 



constant and the factor i is required by the hermiticity of the density matrix. 
Since, however, K ■ q = {p\ this addition does not affect (in the 

momentum representation) the density matrix elements we are interested in 
0. It does affect the deduced size of the interaction region, but this problem 
is not discussed in the present paper. 

Let us make some remarks about normalization. It is often convenient to 
use the normalization 

" d'ppf\p,p) = l. (17) 
In other cases an invariant normalization may be preferable 



-pZMp) = i. (18) 

where Ep = y/m'^ + p^. The relation between the two density matrices is 



ptlMp') = ^/E^'pf\p,p') (19) 

and can be used at any time to go from one normalization to the other. For 
definiteness we will use the matrix pf'^ normalized according to (|l 



The calculation of multiparticle effects can be greatly simplified, when 
for the input density matrix pf*^ its eigenfunctions "ipnip) and its eigenvalues 
A„ are known. These are defined by the equation 

d'p'pf\p,p')Mp') = XnMp), (20) 

The eigenfunctions and the eigenvalues are known both for the pure state 
model (p!3| ) and for the Gaussian model (|I5|). For the pure state model 
obviously 

i^o{p) = i'iP) and A„ = 6n,o- (21) 

with the remaining eigenfunctions constrained only by the condition that 
they are orthogonal to ipo and to each other. For the Gaussian model in one 



dimension 0, [T^ 
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An 

where n = 0, 1, . . ., 



a 



(1 - z)z\ 



2RA - 1 



exp ^ H^{ap), (22) 

(23) 

In the three-dimensional case the index n becomes the set {nx,ny,ny}. The 
eigenfunctions and the eigenvalues are 

V'n(p) = i^nSPx)'^ny{Py)'^nSPz)] K = A„^A„^An,. (25) 

In terms of its eigenfunctions and eigenvalues the input density matrix is 



Pi (p,p') = 1] A„'0„(p)^;(p')- 

n 

The normalization condition implies that 



(26) 



(27) 



4 Step 2 — independent particles stage 

The density matrix p^^^ is also an auxiliary construct and does not correspond 
to an existing physical system. It describes a system of n independent, 
distinguishable particles. The particles are independent in the sense that the 
average of the product of any two single particle operators acting on different 
particles is equal to the corresponding product of averages 



(Ol(Pl)02(p2))n = (Oi(pi))„(Ol(pi))„ 

where the averages are defined by the standard formula 



(28) 



(29) 



One could ask what is the relation between the physical assumption that 
pions are emitted independently and the factorization of the weight function 
in the definition of the density matrix. E.g. in the wave packet picture does 
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p„(pi, . . . ,p„) = |ai, • • • ,an)p„(ai, • • • ,an)(ai, • • • ,an|, (30) 



ai,...,a„ 
n 



p„(ai, . . . = nPi(«fc)' (31) 

k=l 

imply that the emission is independent. As easily checked from definition 
p8| ) the answer is affirmative, if the wave functions (pi, . . . . . . , a„) 

are products of single particle wave functions (pfclo^fc)- This is the case in 
step two of the present model. If the n-particle wave function corresponds 
to correlated particles, the factorizability of the weight function does not 
help. This is the case after the wave functions are symmetrized, since then 
the Bose-Einstein correlations appear. It is a matter of taste, whether these 
correlations are ascribed to the symmetrization of the states 
which are part of the density operator, or to the symmetrization to the 
external states . . . ,p„). The first choice was made e.g. by Zimanyi and 



Csorgo 1 10 1 and the second e.g. in the classical paper of the Goldhabers Lee 



and Pais [|14| . This is not a physical distinction, however, because the scalar 



product defining the n particle wave function does not depend on whether 
the first factor, the second factor, or both factors got symmetrized. 

Thus, the model is an independent emission model [|, ||, ^ in a well- 
defined though somewhat formal sense. When pf^ depends on unobservable 
parameters, independence holds for fixed values of these parameters and 
would be destroyed, if one averaged over them. 

One could define for distinguishable particles, in analogy to step 4, a 
density matrix, which involves all the multiplicities 

P^°^ = i:^e->(°). (32) 

n=0 

It is easy to check that also this matrix corresponds to independent particle 
production. A different choice of the probabilities p*^''^ (n) would correspond 
to correlated production even at the stage when the particles are consid- 
ered distinguishable 0. Of course, introducing correlations by modifying 
the input multiplicity distribution only, without modifying the independent 
production assumption for each given multiplicity, is not the most general 
way of introducing correlations. 
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5 Step 3 — exclusive momentum distributions 



The diagonal elements of the symmetrized n-particle density matrix obtained 
in step 3 yield the momentum distribution for n identical particles, when no 
more particles of this kind have been produced. There are no constraints on 
the production of particles of other kinds. Thus, strictly speaking, this is a 
semiinclusive distribution. Since, however, in this paper particles of other 
kinds are ignored and could just as well be assumed to be absent, we have 
called this distribution exclusive. It is in general not normalized, even when 
the single particle density matrix pf^ is normalized. We will use the notation 
g: Wo = 1 and for ri > 

1 /■ " 

Wn = Trpn = —J2 II Pi\p<jik,PTk)d^Pk (33) 

a,T ■' k=l 

Using the matrix pf^^i instead of pf^ one can make this formula invariant. 
The same remark applies to our further formulae and we will not repeat it. 
For further use it is convenient to define the generating functional [|, H: 
Wq[u] = 1 and for n > 

/n 
Pn{pi, ■■■,Pn) n u{pk)d'^Pk- (34) 
k=l 

For u{p) = 1 one recovers the parameters Wn- Wn = Wn[l]. 

Since each permutation can be decomposed into cycles, the parameters 
Wn for n > can be expressed in terms of the simpler parameters defined 
for A; > by 

/n 
p'f\puP2)p?\p2,P3) . ..p?\pk,Pi) n d'Pr (35) 

If the input matrix pf^ is normalized by (p!7|), Ci = 1. The corresponding 
functional is 

» n 

Ck[u] = I Pi\pi,P2)pi\p2,Pz) ■ ■■Pi\pk,Pi) n '^ipM^Pj- (36) 
with Cfc[l] = Ck. 
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The functional VF„[m] for n > can be expressed in terms of the func- 
tional Ck[u] according to the formula 

WM-r, E n«!. (37) 

ni,...,n„k=l 

The sum is over all the sets of nonnegative integers {ni, . . . ,nn} satisfying 
the relation J2k=i ^^fc = n, or equivalently over all the decompositions of the 
set of permutations of n objects into cycles, so that there are cycles of 
length k. Note that Trp^'^^ = is equal to the term rik = n5k,i of this sum. 
Thus, all the further terms in can be interpreted as corrections due to 

symmet r iz at ion . 

At this point one can write down the normalized n-particle exclusive 
momentum distribution 

P{pi,...Jn) = ^Pn{.Pl,---Pn)- (38) 

The distribution is here normalized to unity, but a change of this convention 
would be trivial. 

Using representation (pB]) for the input single particle density matrix, one 
can write the parameters in the form 

Ck = Y.>t- (39) 

n 

For the pure state model model (|T^ 

Cu[u] = (1 rf='p|^(p)|\(p)'; A; = 1,2,..., (40) 

which implies for each k: = 1 and VF^M = k\Ck[u\. 

For the Gaussian model in three dimensions, after summing three geo- 
metrical progressions, 

C.= \{ (41) 

i=x,y,z i 

Two limits are here of interest |^, |10[. When the phase space per particle 
is minimal, for i = x,y,z: 2Aj_Rj — 1 and Zi 0. Consequently, Ck ^ I 
for all k > 0. The state becomes pure as it should, when due to Einstein's 
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condensation practically all the particle are in the single particle state cor- 
responding to the eigenvalues Aoi- When the phase space is large, for each i: 
AiRi — > cxo and — > 1. Consequently, Ck — > 5fc,i. In this limit multiparticle 
effects become negligible. 

Let us note that in order to calculate momentum distributions one uses 
only the diagonal elements of the density matrix. Therefore, for this calcula- 
tion any density matrix can be replaced by a diagonal matrix with the same 
diagonal elements. For instance, the density operator proposed in refs P, |[ : 

p = e~" exp 

where annihilation and creation operators for particles with mo- 

mentum p, yields a density matrix nondiagonal in n. All that matters for the 
calculation of momentum distributions, however, are the diagonal elements 
and in our notation one finds for each subspace of n-particle states 

p^'\n){p,, . . . ,Pn\Pn\pi, . . . = — n \JiPk)\'- (43) 

k=l 

Thus choosing p^^\n) = exp(— n)n"/(n!)^ one can replace the density oper- 
ator (^) by the much simpler density operator of the pure state case. For 
given multiplicity the probability distribution for the particle momenta is a 
product of single particle momentum distributions. Therefore, particle mo- 
menta are uncorrelated. In particular, no Bose- Einstein correlations are seen 
in the distributions of relative momenta. This corresponds to the coherent 
case of the GKW model If, on the other hand, the unobservable param- 
eters have not been averaged out, the (unnormalized) exclusive momentum 
distribution is 

{p^'\n)pX{Pi, ■ ■ • ,Pn)) = ^e-^ n \-JiPk)\'^ , (44) 

where (...) denotes averaging over the unobservable parameters. This dis- 
tribution is not a product of single particle distributions any more and con- 
sequently correlations among the particles are present. It corresponds to the 
chaotic case of the GKW model. 



dp J (p) at 



|0)(0|exp 



—i J dp J* {p)ap 



(42) 
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6 Step 4 — multiplicity distribution 



Let us note first that, since in general Wn = Trpn ^ 1, the coefficients p^^\n) 
are not the probabilities for producing n particles. These (unnormalized) 
probabilities are 

p{n)=p^^^Wn. (45) 

Substituting the Poisson formula for p^^\n) and formula ( pT] ) with u = 1 for 
Trpn one finds 

'^p{n) = exp 

n 

For non-Poissonian input multiplicity distributions the corresponding sum- 
mation usually must be done numerically In the following we limit our 
discussion to the Poissonian case. The probability distribution pin) can be 
normalized, if and only if the series in the exponent in (^Bf ) converges. This 
is the case if 

V lim l(Ck = z/Ao < 1, (47) 

where Aq is the largest eigenvalue of the input single particle density matrix. 
The series diverges when z/ > Ag^. This happens for the pure state model 
when 



oo 



(46) 



u > 1 (48) 

and for the Gaussian model, when 

J^>^o= n d + ^iRi)- (49) 

i=x,y,z 

In the terminology of thermodynamics u = uq corresponds to a singularity 
of the (grand)partition function, i.e. (usually, see below) to a phase tran- 
sition. We show further that this is Einstein's condensation. Note that we 
have given here a complete derivation of the formula for uq in the Gaussian 
case. Without using the representation of the input density matrix in terms 
of its eigenvalues, this calculation would have taken many pages ([110] and 



references given there). As Zimanyi and Csorgo put it [|T0[ "these solutions 
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are not easily obtained" — and they worked for simplicity with the spheri- 
cally symmetric case, where Aj and Ri do not depend on i. For u < uq the 
normalized probability distribution for finding n particles is 



Pi^) = —rWnexp 
n\ 



oo 



k=l 



(50) 



When u tends to z/q from below, the sum in the exponent of the normalizing 
factor tends to infinity and the probability of finding n < hq particles for 
any finite Uq tends to zero. This effect disappears, when the eigenvalue Aq 
and the corresponding eigenstate ipQ are removed from the density matrix 
([26|). Thus, the surplus of particles condenses in the state ipQ - Einstein's 



condensation occurs. Since in the pure state model all the particles are in 
one pure state anyway, there is no Einstein condensation in this model for 
1. 



For comparison with experiment one needs the multiplicity distribution, 
its moments, cumulants etc, as well as inclusive single particle, two particle 
etc. distributions. All this information is conveniently summarized in the 
generating functional^ 



= Y.p^^\n)Wn[u] = exp 



n=0 



> — —u - V 



.k=l 



(51) 



This formula differs from the corresponding formula in ref. P by a momen- 
tum independent factor, which does not affect the logarithmic derivatives of 
interest. For the pure state model 



(52) 



1 — u J d^p\ip{p)\'^u{p) 

In order to calculate the moments and/or cumulants of the multiplicity 
distribution, one chooses the function u constant. Let us denote this constant 
by z. Then 



Wn[z] = z^W^, (53) 

Ck[z] = z'Ck (54) 

^Many of the steps below become obvious, when one translates the argument into the 
language of thermodynamics |p^ . 
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Let us calculate for example 



("> = (^™) =E-'C.. (55) 

= 1 fc = l 




{n')-{nr = ( ^ ) +(n) = f:z.'=fcC, = z.^(^). (56) 

z = l fc=l 

These formulae can be simplified by using the eigenvalue expansion. For 
instance, assuming that the geometrical series converge, one finds 

(") = Y. (57) 

n=0 J- ~ 

For given p^i \ i.e. for given eigenvalues A^, the parameter u should be chosen 
so as to reproduce the observed value of (n) . With increasing (n) the param- 



eter u increases. For {n) — »• 00, uXq — >■ 1. Each term in the sum (|57|), except 
for the first one, tends to a finite limit. Let us assume that the sum of these 
limits is non zero and finite. Then, with (n) increasing to large values, only 
the first term keeps growing significantly. The k-th term in the sum gives the 
average population of the k-th eigenstate of the density operator. Thus for 
(n) large and growing, almost all the additional particles due to the increase 
of u land in one state {\ipo)) — Einstein's condensation occurs. The formula 
for the dispersion of the multiplicity distribution can be rewritten as 

(»') - = E 7Y^2 (88) 

n=0 U ~ ^-^n) 

In the pure state model, where each of the coefficients Ck equals one, the 
series for (n) converges only if z/ < 1 and then 

{n) = for z/ < 1. (59) 

For this model symmetrization from step 3 reduces to a multiplication of the 
density matrix by n\. The unnormalized probability for producing n particles 

p{n) = e'V", (60) 
which explains the singularity in (n) for u 1. 
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As seen from these examples, when the eigenvalues A„ are known, the 
numerical evaluation of the moments of the particle multiplicity distribution 
is easy. 

7 Step 4 — inclusive distribution of momenta 

When calculating the inclusive distribution of momenta it is convenient to 
use the functions 



k=l 



(0) 

Pi 



ip,p'), 



(61) 



where 



(0) 

Pi 



{P,P') = J d^Pl ■ ■■d^PnPl\p,Pl)Pl'\pi,P2) ■ ■■Pl\pk,p'), 

or equivalently 



(0) 



Pi 



Using the latter notation 



(p,P') = T.MP)^'nrn{p')- 



L{P,P') = J2'^n{p)llj*n{p')z \ 

1 — VA. 



As seen from this formula 



L*{p,p') = L{p',p). 



(62) 



(63) 



(64) 



(65) 



Actually, L{p,p') is the symmetrized single particle inclusive density matrix 
(see below). 

The various inclusive distributions of momenta can be evaluated as func- 
tional derivatives of the generating functional ^[u] aiu — 1. Thus, the single 
particle distribution and the two particle correlation function are 



Pi{p) = (^^ln$Mj =L(p,p), (66) 
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= \LiPi,P2)\'- (67) 

For the pure state model L{p,p') = Pi'\p,p') and there are no sym- 
metrization effects in the inclusive distribution of momenta. For the Gaus- 
sian model, using the formulae from Section 3, one can obtain the function 
L{p,p') as a power series in u. The coefficients of this series are explicitly 
known Gaussians. 



8 Conclusions 

Let us summarize our conclusions and add a few comments. 

• Many models used for the description of Bose-Einstein correlations in 
high energy mutiple-particle-production processes belong to the class 
of factorizable models described in the present paper. As example of 
a model which does not, let us quote the string model developed by 
Andersson and collaborators |]I6| , p!7| , |T8[ . This model is based on a very 
different picture of the particle production process, though it has been 
suggested [|TU] that numerically its predictions might be very close to 



the predictions of a suitably chosen factorizable model. In practice the 
models are often supplemented with corrections for resonance produc- 
tion, final state strong interactions. Coulomb interactions etc. We have 
not discussed these problems here. 

The input consists of an input multiplicity distribution (usually Pois- 
sonian) and an input single particle density matrix in the momentum 
representation. There is a great variety of physical pictures used to 
suggests these inputs, but there is no consensus on which is the best. 

Once the input has been chosen, it is standard to derive general formu- 
lae for the inclusive and exclusive distributions of momenta, multiplic- 
ity distributions and various correlation functions and coefficients. In 
practice, however, it is usually difficult to use these formulae beyond 
the few particle distributions uncorrected for multiple particle effects. 

The multiparticle effects can be more easily included, when the eigen- 
functions and eigenvalues of the input single particle particle density 
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matrix are known. Examples are the two cases, where exact solutions 
are known, i.e. the Gaussian model and the pure state model. In gen- 
eral, Einstein's condensation occurs, when the model is not the pure 
state model and the parameter v of the input Poissonian multiplicity 
distribution is equal to the inverse of the largest eigenvalue of the input 
single particle density matrix. 
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